Competition between zero bias anomaly and proximity effect in disordered systems 
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We investigate the suppression of the proximity effect in mesoscopic normal metal/superconductor 
systems induced by Coulomb interaction. We identify and elucidate the mechanism by which 
disorder leads to an amplification of this effect. In particular, for strong enough disorder, the 
proximity effect is shown to vanish. An expression for the tunneling density of states is obtained, 
and experimental applications are discussed. 
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The influence of a superconductor (S) brought into 
contact with a normal (N) disordered metal has drawn 
great interest recently both theoretically |l|-0| and ex- 
perimentally S . In such a geometry, the superconductor 
tends to impart aspects of its superconducting character 
on to the normal region. Underlying this effect is the 
mechanism of Andreev scattering by which an electron 
from the normal region is reflected at the SN-interface as 
a hole, and a Cooper pair is added to the superconduct- 
ing condensate. This phenomenon has a striking effect on 
spectral and transport properties, known collectively as 
the proximity effect (PE) 0. In particular, if the normal 
region is finite, a minigap Eg develops in the excitation 
spectrum, the size of which depends on the typical time 
between successive Andreev reflections for quasiparticles 
in the normal metal, which in turn depends on the trans- 
parency T of the SN-interface GUsl . 

The minigap in the excitation spectrum can be mea- 
sured through the tunneling conductance of an external 
probe weakly coupled to the normal system. In the ab- 
sence of a superconductor, the repulsive electron-electron 
(e-e) interactions in the normal metal give rise to a sup- 
pression of the tunneling density of states (TDoS) at 
small voltages [^-11 . This phenomenon, known as the 
zero bias anomaly (ZBA), can be understood as the par- 
tial blocking of electron tunneling due to Coulomb re- 
pulsion from its own charge cloud |11| . Since in low- 
dimensional dirty metals the electron charge spreads 
more slowly than in clean metals, this effect is enhanced. 
One therefore expects that it will play an important role 
in small dirty SN structures, where the ZBA is combined 
with the gap induced by the PE. 

Although various important aspects of the role of 
Coulomb interactions in mesoscopic SN (and S) struc- 
tures have already been investigated ra,Qq;na,0[, the 
disorder amplified interplay between 'Coulomb blocking' 
and PE described above has so far not been discussed. 

To be specific, it will be shown in this Letter how the 
mechanism of Coulomb blocking renormalizes the trans- 
parency of the tunnel barrier at the SN interface, thereby 
leading to a suppression of the minigap Eg. In particu- 
lar, we find that, above a critical strength of disorder. 



the minigap vanishes altogether. In addition to affecting 
the size of the gap, the same mechanism also leads to a 
(non-singular) suppression of the TDoS. (Note that, un- 
like in the absence of Coulomb interaction, the position 
of the gap and the enhancement of the TDoS are not con- 
strained by sum-rules.) Although in this Letter the focus 
will be on the behavior of the TDoS, it is worth mention- 
ing that other properties (e.g. the Josephson current) are 
also affected by the suppression mechanism above. 

Such behavior has important experimental conse- 
quences: If the normal metal is relatively clean, the 
Coulomb blocking effect is weak. At temperatures T ^ 
Eg, a measurement of the differential conductance dl/dV 
as a function of voltage V = e/e for a current injected 
from an external probe into the normal part of the SN 
structure reveals an enhancement of the TDoS above a 
sharp minigap edge. As the effective WM strength of the 
disorder in the normal region is increased both the mini- 
gap Eg and the enhancement near Eg are diminished and 
a superimposed signature of the standard ZBA suppres- 
sion in the TDoS appears. At a critical value of disorder 
the PE disappears altogether, and only the conventional 
ZBA remains. 

To quantitatively study the PE in an interacting en- 
vironment, we employ the non- linear cr-model for disor- 
dered superconductors |15|Jl6|| . The conventional quasi- 
classical picture 113,0 of the PE is obtained by evaluat- 
ing the (non-interacting) cr-model functional in a saddle- 
point (SP) approximation |16|. We then combine this 
method with the standard replica field theoretic approach 
to study Coulomb interaction effects in disordered con- 
ductors [Q . Following this strategy, rather than conven- 
tional diagrammatic perturbation theory, is motivated by 
the fact that the latter is inapplicable in PE dominated 
environments (see below). 

Although the theory is not restricted to a particular 
geometry, to be concrete we consider a diffusive normal 
metallic film of width d separated from a clean bulk su- 
perconductor by a tunnel barrier with transparency T. 
In the dirty limit, (e. A) < r^^ <C Ep, the impurity aver- 
aged (replicated) quantum partition function is obtained 
from a functional field integral involving the low energy 



effective action [ |16| , |19|1 

S[Q] = ^ / dr tr{i?(VQ)2 - A{V + S)Q- 

-27<5(z)[Q(z+)-Q(z-)]'}+ (1) 

+ ^fdr dv' tr [V{v)T^\v - r')F(r')] , 

where D denotes the bare diffusion constant, E(r) — 
(A(r),0,e) • <tP'\ e represents the diagonal matrix of 
Fermion Matsubara frequencies e„ = 7rT(2n + 1), t is 
the mean free time, Ep is the Fermi energy, and the 
order parameter A, specified in an appropriate gauge, 
takes a constant non-zero value in the superconductor 
and vanishes in the normal region. In our slab geome- 
try the film is in the xy plane the interface is located at 
z = and the metal at z < 0. The functional integration 

extends over matrix fields Q = T ( A ® trf 1 T~^, where 

(jP'^ denote Pauli matrices operating in the particle/hole 
sector, A — sgn(£), and the generators T have degrees of 
freedom corresponding to the spin and particle/hole sec- 
tors, as well as replica indices and Matsubara frequencies 
| [19[ . The strength of coupling across the SN-interface is 
controlled by 7 = cvpT jiX — T), where T denotes the 
transparency of the tunnel barrier and c represents a nu- 
merical constant of order unity Ig^l- In the following 
we assume that 7/d <^ Ed <^ ^^ where E^ = Djd} is 
the inverse transport time across the normal region. As 
charge transport across the SN-interface is mediated by 
the coupling term c>c 7 in action (|l|), the fields are sub- 
ject to the boundary condition on the normal derivative 
dzQ\z=n = at the SN-interface. Finally, the influence 
of the Coulomb interaction enters the effective action 
through the Hubbard-Stratonovich field V . Here atten- 
tion is limited to the density-density channel which gives 
rise to singular modifications of the SN-coupling. The re- 
pulsive interaction in the Cooper channel suppresses the 
minigap however not in a singular manner [[7|. In deriv- 
ing the low energy effective action, several mechanisms 
of renormalization have already been taken into account: 
Firstly, the repulsive e-e interaction in the bulk super- 
conductor leads to high energy renormalization of the 
implicit electron-phonon coupling, which in turn modifies 
the superconducting order parameter [gl|. Secondly, the 
kernel Fq ^ includes screening effects both from standard 
random phase approximation corrections pO] , as well as 
from nearby metallic (or superconducting) objects p^ . 
Due to the presence of the order parameter A, stan- 
dard perturbative approaches Jlfl] for evaluating the ac- 
tion dl^) fail. Following Ref. [|16|, we therefore begin by 
subjecting the action to a mean-field analysis. A varia- 
tion of the action (|l|) with respect to Q subject to the 
constraint Q^ = 1, generates the SP equation 



DViQVQ) + [Q, S + V] 



j[QiO+)6iz-0~),Q] 
+ 7[Q(0-)(5(z-0+),Q]. 



In the non-interacting theory, this equation is the con- 
ventional Usadel equation |Q describing the spatial de- 
pendence of the average quasi-classical Green function. 
The right hand side of Eq. (||), together with dzQ\z=o re- 
produces the boundary conditions to the Usadel equation 
in the presence of a tunneling barrier between N and S 
pO| . These boundary conditions implies that in the low 
energy regime, e <^ Ed, Q does not fluctuate across the 
N layer. (The minimum energetic cost associated with 
a fluctuating mode compatible with the boundary con- 
ditions is of 0{Ed).) As a consequence the boundary 
i5- function 6{z — 0^) can be replaced by d~^ which sub- 
stantially simplifies the solution of (g). 

We begin by investigating the non-interacting system: 
V — 0. It is straightforward to check that in the low 
energy regime Eq. (g) is solved by Q{z < 0) ~ cr^ ® 1 
H and. 



Q{z > 0) = A..^ = cos 6* (g) C7^ + sine* 



ph. 



.Ph 



(3) 



where the diagonal matrix 9 has elements 6'„ = 
arccos(e„/S„) with E^ = el + {E°)^ and S° = 7/d. This 
result implies the existence of a minigap in the excitation 
spectrum of the normal region of size E2 — 7/d ^ i?d. 
Applied to the TDoS, v = foRe(tr(7'^ ® Is ® erf*" (g) 
A)Q)q/4, where vq denotes the DoS of the bulk normal 
metal, 1^ is the identity matrix in the spin sector and 
V^ a. projector in replica space jlffl, the mean- field solu- 
tion implies a square root singularity at E^. (In passing 
we note that, due to the smallness of the SN-coupling, 
invoking a self-consistent renormalization scheme for the 
order parameter near the interface is unnecessary |23J.) 

Restoring the Coulomb interaction field V, the SP 
equation (0) becomes non-diagonal in Matsubara space. 
An exact solution being now unavailable, a natural ap- 
proach would be to seek a perturbative solution of the 
SP equation. However, as we will see below, such an 
analysis fails to properly account for the renormalization 
of the coupling constant 7. While the former leads to a 
suppression of the TDoS above _£", the latter leads to 
an explicit renormalization of the gap itself. This effect 
parallels the renormalization of E^ encountered in non- 
interacting disordered SN-structures due to mechanisms 
of weak localization (in the particle/hole channel) p6[ . 
Our approach is to include both effects systematically. 

Beginning with a perturbative analysis of the SP equa- 
tion p4[ , we seek a transformation in the normal region, 

Q(0 < z < d) = e-^/2^^gM//2^ {^^ A^} = 0, (4) 

which satisfies Eq. (0) to leading order in both the gen- 
erators W and V. Substituting this ansatz into the SP 
equation and linearizing we obtain 



Dd^W -{W,E} = [V,A^], 



(5) 



(2) 



where the diagonal matrix E = {£"„}. For the slab ge- 
ometry, the solution is given by 



2Vnm 

Wn,n{q) = -^ """, n 2 [(^OS 6'„ - COsOrn) (J3 

t^m + ^n + i->q 



- (sin 



On -sin6',„)cri] ® 1, 



(6) 



where q is the momentum in the xj/-plane and the de- 
nominator has the signifieance of a generahzed diffusion 
pole (which, in the metaUic hmit E^ — > 0, approaches the 
standard form of a diffuson.) 

Substituting Eq. (0) into the action we obtain an ef- 
fective action for the field V , describing how the screen- 
ing properties of the normal metal are modified as a re- 
sult of the diffusive motion in the presence of the su- 
perconductor. In particular, the e-e interaction kernel 
Tq^ -^ r^^-l-II = r~^, where the operator n accounts for 
the anomalous polarizability of the PE influenced metal. 
The dominant contribution to the TDoS derives from 
processes where only small momenta g <C 1/d and ener- 
gies Lo <^ Ed are transferred by the Coulomb interaction. 
In this case screening due to both the superconductor and 
the polarization operator 11 render the Coulomb interac- 
tion effectively short range, r^^(r — r') « A^^, where A 
is a constant. (We remark that a more precise treatment 
of the screening properties may be important in other 
geometries.) 

Using the Coulomb corrected SP solution above to 
compute the zero temperature TDoS, after summation 
over Matsubara frequencies and analytic continuation 
e„ ^ ie + 0+ we obtain the correction psl 



5^{e) 



\Xte{e^ 



E" 



^1/2 



f(x) = ^{X^ - f )-l/2 log [x + (a;2 _ 1)1/2] 



(7) 



where p is the resistivity of the thin normal layer, 
t = {e^ /2TT'^h){p/d) controls the strength of disorder 
and 6{x) = 1 for x > and otherwise. Noting that 
f{x) ~ — log2x for a; ^ 1 and /(I) = 1, we see that 
Eq. (^ reproduces the usual ZBA for e ^ E2^ and sat- 
urates to a finite value at energies close to E^. Signif- 
icantly, no states appear below E'^. From Eq. (0) we 
see that the strength of the square root singularity is 
suppressed by a logarithmic contribution, the coefficient 
inside the log being slightly larger than the coefficient in 
the usual ZBA. There are two effects that distinguish this 
contribution from the ZBA in the absence of the super- 
conductor: Firstly, due to the enhancement of the TDoS 
in the vicinity of _Eg, the single electron propagator de- 
cays (in time) slower than in a normal metal, thereby 
enhancing the Coulomb blocking effect. Secondly, the 
conductance near E^ is larger so that the cloud can relax 
faster, which leads to a smaller interaction correction. As 
can be seen from Eq. (0), the combination of these effects 
does not lead to singular corrections to the TDoS. 

Eq. (0) predicts that the position Eg of the TDoS gap 
is unaffected by interactions, thereby indicating that our 



so-far perturbative analysis is qualitatively incomplete: 
We expect that Coulomb blocking leads to a downward 
renormalization of the coupling strength between N and 
S, and hence to a shrinkage of the gap itself. To quanti- 
tatively describe this effect, we follow Ref. Q, and em- 
ploy a renormalization group (RG) procedure in which 
the fast, high energy degrees of freedom renormalize the 
slow. 

Eq. (M) shows that the main (logarithmically divergent) 
contribution to the TDoS stems from energies larger than 
the (renormalized) minigap Eg where the system behaves 
effectively like a metallic one. This observation dictates 
the strategy of the RG approach: Employing the (renor- 
malized) Eg itself as an infrared RG cutoff we assume 
that for energies e > Eg, the system can be approximated 
by a conventional 2D disordered conductor. For energies 
s ^ Eg the treatment is accurate, while for e ^ Eg we 
expect corrections which are small in the limit Ed ^ Eg 
p6| . In the ultraviolet, the RG procedure is cut off by 
the "high" energy scale Ed- 

The separation of energy scales above allows us to ap- 
ply, without technical modifications, the standard RG 
procedure for 2D disordered systems pC| in the interval 
Eg < {Dq^,e) < Ed- Conceptually, the only difference 
is that the perturbative RG is superimposed on a back- 
ground of the Usadel solution. Carrying out the RG pro- 
cedure, yields a renormalization of the coupling strength 

7, rari 



^(e)=^oe-(*^/4)iog(B,/e)^ 



(8) 



Combined with the relation Eg = ^ /d we obtain a self- 
consistent equation for Eg, 



Eg^E"e-°'^°«'^''/'^<^, 



(9) 



where a = tA/4. The solution of Eq.(y) represents the 
renormalized energy gap in the normal metal layer. 



E* = 



K{if 



a < 1. 



(10) 



Hence, the gap decreases monotonically with increasing 
a (or disorder t) up to a critical value a — 1 above which 
E* vanishes altogether. 

Treating the tunneling from the probe on the same 
footing (i.e. renormalization above Eg as in the case 
of a normal metal, and applying a cut-off at Eg) leads 
to an exponentiation of the first order correction to the 
TDoS |To,ni|. Thus, the asymptotic regimes above and 



below Eg can be summarized as 



1^0 



£g-alog(_Ed/e) 1-^2 





Eg{e) = E^e-°''^°s^^''/'\ 



' Elie) 
Eg{e = 



1-1/2 



E* 



e <e: 

--E*. 



Before concluding, a few additional remarks are in or- 
der. First, in the RG treatment employed here, the ex- 



change and Cooper channels, as well as the renormaliza- 
tion of the diffusion constant and frequency by interac- 
tion and weak localization effects |1^ where neglected. 
Since there is no reason to expect that these terms are 
small in the strong disorder case when \og[Ed/e\ > 1/t, 
the accuracy of the present theory is limited. Since we 
expect all the terms we neglected to further suppress the 
PE, we believe that their effect will be to decrease the 
critical strength of disorder below a — I. 

Secondly, it is instructive to interpret the Coulomb 
renormalization of the proximity effect discussed above 
from an alternative point of view. The fields Q can be 
subjected to a quantum mechanical gauge transformation 
Q -^ e**(5e~**. An interesting reformulation of the the- 
ory arises if $ is chosen so as to eliminate the most singu- 
lar momentum/frequency components of the disorder en- 
hanced Coulomb potential V. As for the case of extended 
N-systems, it has recently been shown g^'lll ^^^^ gaug- 
ing the theory in this way, leads to the nonperturbative 
expression for the ZBA ||l^. In the problem considered 
here, the situation is different, the gauge transformation 
not only removes the singular components of the inter- 
action but also directly couples to the phase of the or- 
der parameter in the superconductor. More specifically, 
gauging the theory as outlined above, induces dynamical 
boundary phase fluctuations. These effectively reduce 
the SN-coupling and hence suppress the proximity effect. 
These considerations indicate that the gauge transforma- 
tion p^,p8[ effectively mediates between the phase and 
the charge representation [E9| of charging phenomena in 
disordered SN-j unctions. 

Finally, the suppression of the PE due to e-e interac- 
tions is expected to be present in other hybrid SN struc- 
tures as well. In particular, corrections due to Coulomb 
blocking would be expected to be larger when the dimen- 
sion of the normal metal is smaller than two. We expect 
that the main features of the analysis remain valid: the 
suppression of the PE with increasing disorder, and its 
disappearance at a critical value. 

Summarizing, we have studied the interplay between 
the PE and e-e interactions in disordered SN-systems. 
Focusing on an analysis of the TDoS, it was shown how 
disorder amplifies the destructive effect of interactions on 
the PE. Above a certain disorder concentration, the sup- 
pression is complete in the sense that the characteristic 
PE-induced minigap structure in the TDoS disappears. 
It was argued that the phenomenon should be of concrete 
relevance for the tunnel-spectroscopic observation of the 
PE in disordered mesoscopic SN junctions. 
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